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Abstract 

We study Galileon theories that emerge in ghost-free massive gravity. In 
particular, we focus on a sub-class of these theories where the Galileons can 
be completely decoupled from the tensor Lagrangian. These Galileons dif- 
fer from generic ones - they have interrelated coefficients of the cubic and 
quartic terms, and most importantly, a non-standard coupling to external 
stress-tensors, governed by the same coefficient. We show that this theory 
has no static stable spherically symmetric solutions that would interpolate 
from the Vainshtein region to flat space; these two regions cannot be smoothly 
matched for the sign of the coefficient for which fluctuations are stable. In- 
stead, for this sign choice, a solution in the Vainshtein domain is matched 
onto a cosmological background. Small fluctuations above this solution are 
stable, and sub-luminal. We discuss observational constraints on this theory, 
within the quantum eff^ective Lagrangian approach, and argue that having a 
graviton mass of the order of the present-day Hubble parameter, is consistent 
with the data. Last but not least, we also present a general class of cosmologi- 
cal solutions in this theory, some of which exhibit the de-mixing phenomenon, 
previously found for the self-accelerated solution. 



1 Introduction and summary 



In the present work we will be interested in those Galileons that describe helicity-0 
modes of tensor fields: In the context of the DGP model pj the cubic Galileon was 
found in Ref. p], and was shown to describe properties of the helicity-0 mode. More 
general, quartic and quintic Galileons, were studied in [3j. 

Most relevant to the present work, is the fact that all the Galileons appear in 
the Lagrangian for the helicity-0 mode of a massive graviton [U |5] (see j6] for earlier 
evidence). Although the massive gravity Galileons are similar to the generic ones 
[3], there is a nuance that distinguishes them: the massive gravity Galileons come 
with specific coefficients, and most importantly, with a novel coupling to the matter 
stress-tensor governed by the same coefficient [SI IH E]- 

One motivation for this work is to study spherically symmetric solutions in this 
special class of massive gravity Galileons. In particular, in the present work we will 
focus on the sub-class of these theories for which the helicity-0 and helicity-2 modes 
can be diagonalized [4j. We refer to this sub-class as Restricted Galileons. 

Spherically symmetric solutions in general Galileon theories are well understood 
[2]: generically, one gets a solution that transitions from the Vainshtein region to 
the asymptotically fiat one; furthermore, typically, stable fiuctuations, both in and 
outside of the Vainshtein region |8], exhibit super-luminal propagation [3]. This 
however, does not necessarily lead to acausality - attempts to create closed time- 
like curves bring one beyond the domain of validity of an effective field theory [D]. 
Pending a classical, or quantum, completion of these theories at the strong scale, it 
makes sense to postpone the discussion of this issue. For a review of other interesting 
theoretical constructions with Galileons, see Refs. [10]. 

We will show that, there are significant differences in the case of restricted 
Galileons (ReG). In particular, we will prove that in ReG no stable static spherically 
symmetric solutions exist that would interpolate between the Vainshtein region to 
an asymptotically flat domain. Instead, the Vainshtein region is naturally matched 
onto a solution that does not asymptote to flat space, but does to a cosmological 
background. The perturbations on this solution are stable, and sub-luminal! 

For the exact static solution, we address the issue of the bound on the graviton 
mass, as inferred from it. We show that having the graviton mass as small as the 
value of the Hubble parameter today is consistent with the bounds derived from the 
table-top sub-millimeter measurements of gravity-competing forces. Our conclusions 
differ from those in [11], and we discuss the origin of this difference. 

The exact static solution, if applied to distribution of sources like the ones in our 
Universe (galaxies, clusters, etc), is consistent with the observations, due to the fact 
that the Vainshtein radius in such a universe gets pushed - statistically speaking - 
closer to the Hubble scale. There may however be isolated clusters in the universe 
for which beyond-the-Vainshtein-scale physics can be explored. 

One could ask what would happen to a dilute distribution of matter that initially 
has no Vainshtein region, which is placed in an asymptotically flat space, and is 
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being adiabatically collapsed to form a smaller object with higher densities. Then, 
the energy density and pressure in the nonlinear helicity-0 field blow up when the size 
of the distribution becomes of the order of the Vainshtein radius, preventing such an 
adiabatic collapse. Thus, we'd expect the respective spherically symmetric collapsing 
solution, if it exists at all, to have time dependence that is fast as compared to the 
Vainshtein scale. In Appendix A we also present a class of more general time- 
dependent solutions, however, non of them reach an asymptotically fiat space for 
any constant time slice. 

Three important comments on the literature: 

(A) All the results of the present work are obtained in the decoupling limit. As 
to whether our results apply to the solutions of the full theory, discussed in Refs. 
DSl [131 [II ESI CBj, remains to be seen. 

(B) Refs. |17l [T8] have recently shown the existence of superluminal shocks 
for some parameter space of massive gravity. Such superluminal solutions, like in 
Galileon theories, appear well-within the effective field theory domain. However, as 
briefiy discussed above in the 4th paragraph of this section, and in detail in Ref. 

the question whether the corresponding solutions can lead to acausal stable 
closed timelike curves, needs to be addressed including perturbations above these 
solutions, in the context of an effective field theory. In cases considered so far, 
the fiuctuations above such solutions bring one beyond the effective field theory [9| 

- one more argument that these theories need either calculational prescription, or 
classical/quantum completion below/at the strong sea 

(C) Ref. [I?] has also brought up 1970 lecture notes by Zumino [19], where Wess 
and Zumino had proposed four models of a massive spin-2 hadron, /^j,, interacting 
with Einstein's gravity for g^^. Using the recently developed methods in [5| 120 } [2T | 
[221 [23], we can tell now that the two (interrelated) "f-g" models of Wess and Zumino 
are ghost-free as full non-linear theories and represent a subset of more general 
ghost-free bi-gravities if reduced to a single massive graviton, they represent a 
subset of the ghost-free massive gravities [S], from which [21] originates. The two 
other theories proposed in [19] have Boulware-Deser (BD) ghosts [23]. None of the 
conclusions on the absence or presence of the BD ghost can be deduced from Ref. 
[19j itself, not surprisingly, as it precedes the BD work [24|. More details will be 
given elsewhere. 

^One example, where super luminality was claimed to be removed, at least for spherically sym- 
metric sources at the classical level, is a bi-galileon theory |2l] . This requires a low energy remnant 

- the second Galileon which is massless and operates at arbitrarily low scales. We note that certain 
bi-galileons naturally emerge in generalization of ghost-free massive gravity |1] to a theory with a 
dilaton-like particle, quasidilaton |^. Although the issue of superluminality for the quasidilaton 
has not been studied yet, predictions of this theory are not brought in conflict with observations 
by the presence of the low-energy remnant, the quasidilaton |26| . 
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2 The Theory 



In the decoupling limit of massive gravity, the Lagrangian describing helicity-2 
states, contained in h^^, and helicity-0 state, denoted by vr, with their coupling 
to matter stress-tensor T^^, is given by 

^ = -Ih^'S^lu'h., + h^'' (x« + ^Xj^J + ^Xj^^ + -^h^'T,., (1) 

where, the three conserved symmetric tensors xj^i}{Il) depend on the second deriva- 
tives of helicity-0 field II^i, = d^dyTT, and can be cast as follows [27j . 
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Ji^(3) ^ 1 ap7p /3<^^rr TT TT 
^^f_Lii 2 -'■-'■op-'--'-pcr-'--'-7T; 

where is the usual Levi-Civita symbol. 

Under the linear diffeomorphisms 5h^j, = d^C,v + d^C,^, the full nonlinear La- 
grangian (II]) is invariant up to a total derivative, while it is exactly invariant under 
the field space galilean transformations, ^vr = v^x'^, with constant . 

Due to the specific structure of this Lagrangian, the coefficients a and /3 do not 
get renormalized by quantum loops; hence, any choice of their values is technically 
natural |29] . 

For generic values of these coefficients the helicity-2 and helicity-0 modes mix, 
as seen from ([1]); this mixing cannot be undone at the full nonlinear level, unless 
/3 = [4J. The choice /3 = corresponds to selecting a special relation, = —a^/A, 
between the coefficients and of the cubic and quartic K terms of the full 
massive theory [Sj. 

For /3 = 0, which as already mentioned is a technically natural choice stable 
under loop corrections, there exists an invertible field redefinition 

— ct 
h^u h^i, + 7rr]f,^ + — gS^vr^j^Tr, (2) 

that decouples the tensor and scalar modes from each-other [4] . The Lagrangian ([1]) 
then takes the form C = C]^ + C^^, where Cji is the linearized Einstein-Hilbert action 
for h minimally coupled to T^^, while 

= I'Dn + lisiO'fo^ + 1 - |n]^) 
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Here, square brackets [.] denote the trace. Note that the absolute value of the only 
free parameter a is immaterial, as it can be absorbed into the scale A3, however, the 
sign of a will be crucial [rescaling A3 would correspond to rescaling of the graviton 
mass; for comparisons with massive gravity, however, we will manifestly keep the 
parameter a, and unrescaled A3 in what follows]. 

Thus, for /9 = 0, the dynamics of the helicity-0 mode is similar to the cubic 
and quartic Galileons, but with some special coefficients in front of these terms; 
moreover, the last term in the second line in ([2]) is a coupling to stress-tensor that is 
not usually considered in the Galileon theories, which however is present in massive 
gravity [6llll[7], (it will plays a crucial role). As mentioned in the previous section, 
we refer to ([3]) as Restricted Galileons, and discuss this theory in the present work. 
The case with /3 7^ 0, which appears to be quantitatively different, will be discussed 
in Ref. [301 . 

As to the parameter a, a 'priori it can have the either sign. However, the last 
term in ([3]), leads to classical renormalization of the kinetic term for vr. For a 
negative sign of a this renormalization is negative, and for most of the reasonable 
sources it would overshoot the sign of the tt kinetic term. Furthermore, we will 
show in Appendix A, that classical renormalization of the vr kinetic term due to the 
non-linear terms in ([2]), although positive, is in fact subdominant. Therefore a < 
case is not physical: when a localized source (with density greater than AgMpi) is 
approached from far away, at some point the vr field will have vanishing kinetic term, 
signaling an infinitely strongly coupled regime; at yet shorted distances the kinetic 
term would fiip its sign (if one could pass into this region), to convert tt into a ghost. 
For this reason, our physical choice will be a > 0. 



2.1 Spherically symmetric solutions 

In this work we study a gravitational field configuration created by a source of finite 
size i?, and a uniform density p. The pressure of the source will also have some 
important consequences, and will be discussed in Appendix B. 

Without the loss of generality, the static spherically symmetric configuration can 
be parametrized by the following ansatz for the metric perturbations 

hm = air), hij = f{r)5ij, (4) 

while for the helicity-0 we begin by assuming the radial ansatz vr = 7r(r). Then, the 
equations of motion for graviton reduces to two ordinary differential equations (this 
is after integrating them once and requiring the solution to vanish at the origin) 

rf' = -l^ + Ayx(i - a\), ra' = - A^X , (5) 

•' Mpir ^ ^ ^ Mpir ^ ^ ^ 

where above, and in what follows we use the notations for A and the Vainshtein 
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radius r*: 



(\ 1/3 
M \ 

and the prime denotes differentiation w.r.t. the radial coordinate r. 

The closed form of the equations of motion for the helicity-0 mode can be ob- 
tained by integrating out a and /, using The net result reads as follows [28j 

3X-6aX' + 2a'X'=l ^{^l . (7) 

I 2fe)' r<R 

We could have also derived this latter equation by simply integrating the Galileon 
equations of ([3]), with the same initial conditions. 

It is now straightforward to show that for a > equation ([7]) has no solution that 
could interpolate from the Vainshtein region to an asymptotically flat one, i.e., to the 
region where A — )■ 0. To see this we look at the cubic polynomial in A at the Ihs of (ITj). 
The necessary and sufficient condition for the solution of ([7]) to interpolate between 
the Vainshtein region and the asymptotically fiat one, is that for the polynomial on 
the Ihs of ([7]) to have only one zero, at A = [3J. However, this condition is not 
satisfied if a > 0: there are three real zeros in this case, Ai = 0, A2,3 = (3 ± \/3)/2a. 
This means that the solution that starts off in the Vainshtein region, where A 3> 1, 
matches onto a solution that at r ^ oo tends to A2 = ■ The latter implies that 
asymptotically n — )■ A3r^A2. 

Then, the solution to ©-([Tj) for a > 0, in the leading approximation, is: 

Outside the Vainshtein radius 



2^ 



A - ^ , a~ |-A, /~^A(l-aA). (8) 



Inside the Vainshtein radius 

r 1 r 9M 

^^^ + - + :w7^' a^/^--— + 0(r2A). (9) 

a^'-^r a la^'-^r^ Mp\r 

This solution describes the static spherically symmetric source of mass M on a 
cosmological background, with the equation of state p/ p = —2/3 + ~ —0.1. 

Next, we would like to proceed with the analysis of fluctuations around this 
solution. As it is clear from the above, we gave up the asymptotic flatness. Therefore, 
it is natural to begin by studying the region far away from the source. Due to the 
decoupling between h and vr we can analyze their fluctuations separately. Those of h 
are linear GR perturbations on a cosmological background with the above-mentioned 
equation of state, taken in the high momentum approximation (momenta being 
larger than the Hubble parameter, which in this case is ~ m). These fluctuations 
are stable, and there is no novelty in this part. 
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Let us then look at the fluctuations of vr, which we decompose as vr = TTci{r) + 
a{x,t), where 7rcz(r) denotes the classical background just found in (|H]). Appendix A 
gives an expression for the Lagrangian of fluctuations about a general background $; 
substituting $ = tTcI into (lA-Iip . we obtain the quadratic Lagrangian for fluctuations 
around (Ml: 



C = 3 



5 + 3^3 



1 + V3 



{d,a) 



(10) 



Hence, the cosmological background is stable against all linear perturbations. More- 
over, the fluctuations propagate with sub-luminal velocities. 

As a next step we study fluctuations in the presence of the source of mass M. 
Outside the Vainshtein radius the just-studied cosmological background dominates, 
and hence the fluctuations are stable and sub-luminal. However, inside the Vain- 
shtein radius the source defines the background and its effects are dominant. The 
fluctuations of h coincide, to a good approximation, with those of linearized GR on 
a background of a static spherically symmetric source. The fluctuations of vr, how- 
ever, need an extra caution. Using again (lA-Iip . we find that in the leading order in 
r/r^, (which is a good approximation in the Vainshtein regime where (r/r^,) -C 1), 
the small excitations around ^ are described by the following Lagrangian: 



£ = 3 



a'" (^)' + 2a"''j\ {d,af - 3a'/' (^)' (5.^)^ - ^(9^a)^ (11) 



Remarkably, radial fluctuations propagate with a slightly sub-luminal speed -the 
deviation from the luminality being suppressed by r/r*. The angular speed of sound, 
on the other hand, is strongly suppressed. This suppression, for a generic quartic 
Galileons, was already emphasized in Ref. [3J. 

It is interesting to ask what happens with a dilute distribution of matter with 
some size R, and density below A|Mpi, placed on an asymptotically flat space. Such 
a source has no Vainshtein region, since r* < R, and there is a static solution with 
decaying vr field at infinity. However, if the source is adiabatically collapsed to form 
a smaller object with density above A|Mpi (i.e., with r^, > R), then the energy 
density and pressure in the nonlinear n field blows up when r^, ~ R, preventing such 
an adiabatic collapse. 



3 Observational prospects ? 

In this section we will consider a hypothetical universe described by the massive 
gravity with /3 = 0. This is a technically natural choice, as it's not ruined, to a 
high degree of accuracy, by loop corrections ^9\. We take the solution (|8]) as our 
background in an empty space, anticipating that localized sources will match onto 
this solution, and not to a solution with flat asymptotics. 
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3.1 Nonperturbative effects in Solar System 

Since the background ([8]) has an unusual equation of state, we will require this fluid 
to be a subdominant component in the universe [§. This can be arranged by choosing 
the parameter a to be of order ten or so (for a flxed graviton mass m ~ ifo)? which 
again, is a technically natural choice j29]. Then, to obtain realistic universe we would 
have to introduce dark energy, and matter density, roughly as we do in the ACDM 
model. Therefore, from the point of view of getting dark energy from modifled 
gravity, the /3 = theory is not interesting, however, our goal here is to understand 
if such a theory, in principle, could be consistent with observational data. 

To this end we would like to discuss observational constraints due to the small 
corrections to the GR results within the Vainshtein radius jBl From the 

solution ^ we easily deduce that the fractional increase of the Newtonian potential 
inside the Vainshtein radius is 

(1 = — - , (12) 

Tglr 

where by we denote the Newtonian potential. This deviation from the standard 
potential will give rise to an additional perihelion advancement of orbiting bodies 
(an additional to the GR effects). For careful account of these effects, within a 
relativistic theory, see [STj [34J. 

Then, the question is whether such deviations can be measured. One of the most 
precisely measured trajectory is that of the Moon orbiting the Earth For the 
Earth-Moon system, and for m ~ iJo; the fraction (fT2|) is of order qEarth-Moon _ 
10~^^/a^/^; this is four orders of magnitude below the same fraction in the DGP 
model [311 l32j (see, also [33]): as such its influence on the perihelion advancement 
of the Lunar orbit is completely negligible and is unlikely to be measurable by the 
near future Lunar Laser Ranging Experiments |35) . 

Next, we discuss whether fluctuations of the vr fleld, above its background, can 
give rise to additional forces competing with gravity. There is an important point, 
in these considerations, as compared to conventional Galileons. The kinetic term 
for the fluctuations a = — iXci gets additional classical renormalization 

where the factor is due to the cubic and quartic Galileon terms evaluated on 
the solution; we have also manifestly shown couplings of the fluctuations to the 
stress-tensor. Note that the kinetic term in (fT3|) gets enhanced by both a factor 
due to the classical background, Z^, and by a source-dependent term aT (both 
of these are huge factors for realistic sources, but as shown in Appendix A, the 



^Although, other more interesting options may also be possible when the density of this fluid 
adds in a meaningful way to the matter density, while its small negative pressure, adds a bit to 
the negative pressure due to dark energy. 
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enhancement due to the aT term dominates within a source). As a result, the 
static force due to exchange of a, between any two objects, say two small metallic 
balls, is suppressed by a tremendous factor determined by the ratio of density and 
pressure in the balls over the density A|Mpi; the latter happens to coincide with the 
critical density, pc — 10'"^^ g/cm^, if m ~ iJo — 10~^^ eV. This fraction is enormous, 
[Pmetaii / -^sMpi) > 10^". As wc wiU discuss in the next section, this suppression 
makes the aT coupling unobservable. Therefore, the exchange due to a single a 
field cannot be gravity-competing force. In the next section we turn to quantum 
effects, and the question whether they could ruin the above conclusion. 

3.2 Comments on quantum effects 

In this section we would like to discuss quantum loop effect. Before we do so, one 
should point out that we're dealing with a theory that is non-renormalizable by 
power-counting, and at least in this sense, requires a completion at the scale A3. 
However, judging from the properties of the theory below A3 this completion, unlike 
in massive non-Abelian vector fields, is not expected to be a conventional one (see 
[5^ I29|. and references therein). On the other hand, discussions in the literature, 
as well as below, are based on an assumption of a certain conventional generic 
completion. This is an assumption that may not be a right one, and therefore, all 
such discussions are on a shaky ground. 

Putting the above concern aside, we would like to analyze the viability of the 
bound imposed on graviton mass in Ref. jTT], based on the conventional effective 
field theory considerations of the decoupling limit Lagrangian To do so we first 
briefly recall what do the quantum corrections lead to in the decoupling limit of 
massive gravity [29j: the unambiguous result, independent of the uncertainties of 
the previous paragraph, is that the terms presented in ([1]) do not get renormalized! 
However, in a conventional quantum effective field theory approach new terms may 
get induced in the IPI action. These new terms include ambiguous power-divergent 
terms, as well as log divergent pieces that carry information about the forward 
scattering of the quanta, and have to be included in the IPI action. What was 
argued for the cubic [37], as well as generic Galileons [3], is that these terms, in the 
conventional approach, and once calculated on a classical background, end up being 
suppressed by a scale of the classical background itself. 

The classical renormalization f lT3|) increases the effective energy/momentum scale 
of strong interactions of the fluctuations in the IPI action [37J. The naive analysis 
shows that the loop-induced terms are suppressed by the effective scale schematically 
written as follows: 



This is valid when the enhancements of the kinetic and gradient terms are of the 
same order; however, (HM fails to describe the strong coupling scale in case of the 




(14) 
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hierarchy among them. For the illustration, consider the following toy model, that 
captures essential features of our theory on a static background: 

C = a\dtTif - b^djirf + ^{dTryUTT, (15) 

where a, b and 6* are assumed to be constants for simplicity. The kinetic term 
for the field tt needs to be canonically re- normalized, vr — )■ vr/a. This results in 
the naive Ag// ~ a6* A3 estimate for the effective cut-off, in the spirit of (IHl). 
However, by calculating the 2 — )■ 2 scattering amplitude one can see that the strong 
coupling scale differs in different channels; this is because the propagator has the 
denominator, — (6/a)^A;^. In fact, for the forward scattering and 6 <C a, the 
M-channel is the strongest. It gives rise to the Ag// ~ (6/6*)^/^a^/^A3 cut-off, while 

— 1/3 

the s-channel amplitude becomes strong only at a6* A3 scale; the latter being 
considerably larger than the former. 

Higher order tree-level diagrams give higher effective scale: the vertices have 
more powers of inverse a than the propagators can compensate for due to the {h/af' 
factor in the diagrams that aren't automatically zerqj. All this remains to be the 
case even after the quartic Galileon is included in ffTSj) . Hence, the effective lowest 
strong coupling scale remains to be: 

Ae// = a^/^(^J A3. (16) 

At this scale, interactions of the fluctuations would become non-perturbative. Let 
us now discuss how this squares with the fact that gravity-competing forces have 
been excluded down to the scales of order 100 microns or so [40j. First, recall 
that in the DGP model, A^j^ is of order a centimeter j3^, the scale much bigger 
then 100 microns. However, as pointed out in |37] this is not likely to be an issue, 
since these fluctuations are very weakly coupled to matter sources, and thus, their 
effects are unlikely to be visible in the submillimeter measurements, even though 
they themselves self-interact strongly at scales below 1cm. 

For the model considered in the present work, it was argued in Ref. |11) that 
the A~yj, is much bigger that a millimeter, if the graviton Compton wavelength is 
taken to be of the Hubble size. Does this mean that such a small graviton mass is 
ruled out by laboratory measurements, as claimed in [llj ? 

To address this question in detail, consider a device measuring putative gravity- 
competing forces - a torque pendulum of the Adelberger's group experiment |40] . 
Then there are two important points: 

(1) Within the matter that the pendulum plates are made of the effective strong 
coupling scale differs significantly from its value outside of the plates. Using (ITB]) . 

^The power divergent terms in the loops, that could induce lower cutoff, are removed by the 
counterterms as per the general philosophy adopted in this approach; the log divergent terms 
however, are related to the tree-level diagrams by unitarity and analyticity. Hence, it suffices to 
get the lowest cutoff of the tree-level diagrams as we did above. 
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and the fact that in this case a ~ 10^^, 6 ~ 6*, we find that in the plates ^^jf — 
10~^° 1000 km ~ 0.1 mm. This is smaller than the thickness of the upper and lower 
pendulum plates, that respectively are about 1.8 and 7.8 mm 

(2) Thus, in the torque pendulum plates, the fiuctuations are weakly coupled 
to themselves, and they can be well- described by the classical Lagrangian down to 
0.1 millimeters. On the other hand, as shown in the previous section, this classical 
Lagrangian gives rise to a tremendous suppression of the coupling aT/M^i, that 
determines how strong/weak the coupling of the fiuctuations to the plates could be. 
For the metallic plates this suppresses the aT/Mpi vertex by an additional factor of 
a-i ^ 10-151 

Therefore, a single a fiuctuation cannot efficiently be emitted by the plates (or 
any realistic source bigger that 0.1 millimeter, for that matter). Even though the 
fiuctuations in the vacuum between the plates become strongly interacting with 
themselves, and could form some bound states of a, what's important is that this 
exotic strongly self-coupled sector cannot couple efficiently to the measuring device 
- the latter coupling due to the exchange of a is a~'^(h/a)~'^ ~ 10-^* of the strength 
of gravity. Therefore, the exchange due to a cannot be gravity-competing force - in 
spite of the fact that the a field self-interacts strongly outside of the plates. This 
suppression was not taken into account in Ref. |11] . making the bound imposed on 
the graviton mass by that work unwarranted. The valid phenomenological bounds 
on the graviton mass are given in Ref. [42j. 

Note that for /3 7^ 0, which was also considered in |llj . the effect of the suppres- 
sion of coupling to matter, due to the d^ad'^ aT^y term, remains valid. However, 
there are technical and conceptual differences for /3 7^ [30], that obfuscate the 
results of Ref. [11] for that case; this and related issues will be presented in detail 
in [30|. 

In conclusion of this section, and to reiterate, we note that many statements 
about the quantum effective theory are based on assumptions about what the UV 
completion might be. It is clear from various standpoints (see, e.g., [29] and refer- 
ences therein), that theories of massive gravity need a completion, or alternatively, 
a dual formulation (perhaps along the lines of [381 EH]), at and above the A3 scale. 
However, judging from the properties of these theories (such as e.g., nonrenormaliz- 
ability of some of the couplings [29j) the putative completion should not be expected 
to be of a conventional and generic type, as preassumed in all the above discussions. 
Therefore, conclusions deduced from such considerations [IT], even if they were plau- 
sible, would not be set in stones. 

"^This scale is greater than the width of the plates, 0.3 microns, used in the experiment of Ref. 
[JT] that exclude forces 14,000 time stronger than gravity at 10 microns. However, the coupling to 
matter that we're discussing is 10""'^® of the coupling of gravity at 0.1mm (see below), and it would 
be unreasonable to expect for it to become 14,000 times stronger than gravity at the micron scale. 
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Appendix A 



In order to study the stability, we split the scalar degree of freedom into the back- 
ground $ and the fluctuation a as follows 



$ + (T. (A-I) 

As a result, to the second order in perturbations, the Lagrangian ([3]) become^ 



+ 



a 



a 



xdnaducr 



^^3 ^^3 
1 



-aT + 



a 



pi 



(A-II) 



For the static clump of dust of constant density, that is for T^j, = p6^6^6{R — r), we 
obtain the following kinetic term (to the leading order) inside the source 



a 



P 



MpiAi 



R 



(A-III) 



Here, the second term in brackets comes from the quartic Galileon term after we have 
substituted the classical solution within Vainshtein region ([9]). After the replacement 
p M/R^, ( lA-IIII) leads to the conclusion that in a < parameter space the scalar 
perturbations have ghost-like kinetic term within the source of radius R <ti a^^^r^,. 
The latter condition can be translated on the language of the source density as 



ap > MpiAg 



(A-IV) 



For phenomenologically interesting value of A3 ^ ~ 1000 km the above bound be- 
comes ap > 10~^^ g/cm^. 



^Here we give only the scalar part of the Lagrangian, since the tensor mode is completely 
decoupled and propagates according to the linearized Einstein-Hilbert action. 
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As we've shown, for a > 0, the asymptotically flat solution cannot be continued 
into the Vainshtein region. However, there are cosmological solutions that are free 
of the above difficulty. 

In Section 2 we have already presented the static solution which asymptotes to 
the cosmological background, rather than the Minkowski space. Here we consider 
the time- dependent ansatz: 

^ = l^ie + Mrl (A-V) 
The introduction of non-zero c modifies ([7]) in the following way 

3(1 + 2ac)\ - 6(a + a^c)\^ + 2a^\^ 

,3 



2(^) (l + 2ac) + c, r>R 
2 (^)^(l + 2ac) + c, r<R 



(A-VI) 



The above-mentioned static cosmological background corresponds to c = in (1 A- VI) . 
However, there is a plethora of available backgrounds, which correspond to different 
values of c and the only discriminating principle among various backgrounds must 
be their stability. In this letter, we do not pursue the detailed analysis of all the 
possible solutions. Merely, we would like to present one more solution which in our 
opinion is quite special. 

From (lA-Vip . it follows that there is a peculiar value of c = — 1/(2q;), which 
makes the equation independent of the mass of the source. As a result, eq.( lA-Vll) 
has three solutions. One of them being A = l/(2a), which corresponds to an unstable 
self-induced de Sitter space found in [27] (it propagates a ghost) and we discard it. 
The other two are quite interesting 

in particular the pion cloud has zero pressure and positive energy. The easiest way 
to see this, is to look at the metric itself 

(A-VIII) 

The effective pressure takes the following form 

P = = -^er'e^'%d^d,h,^ = ^A(/ - a) = 0. (A-IX) 

Moreover, the energy density of the pion fluid, fllling the space, is given by 

p = MpiGoo = -MpiA/ = — > 0, for a > 0. (A-X) 

2a 
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On this background, the quadratic Lagrangian for perturbations is 



(A-XI) 



everywhere in the space. As it is easy to see, the fluctuations are stable and propa- 
gate with the half speed of light. 

We would like to emphasize that the tt profile is independent of the source 



4a 



+ ( 1 + V3 ) r' 



(A-XII) 



This means, that the presence of the static, spherically symmetric and pressure-less 
source does not excite the longitudinal mode of the massive graviton. 



Appendix B 

Here, we would like to discuss the effect of pressure {p ^ p) on the background and 
its stability. The reason is simple, most of the sources in the universe possess small, 
yet nonzero pressure. 

For simplicity, we treat only c = case here, since the generalization is straight- 
forward. The equation of motion for the longitudinal mode ([7]) takes the following 
form 

3A - GaA^ + 2a^\^ 

2 (^) Outside the source 

2 (^)^ 1 + ^ (aA — 3) Inside the source ' 

where Tq is an integration constant, which is determined by the matching condition 
at the surface of the source. As expected, the only change in the field configuration 
is the effective mass. 

Let us recall that we have chosen the sign of a to be positive to avoid ghosts 
inside the source. The danger was coming from the last term of ( ]A-Iip . But, if 



we avoid the ghost by appropriate choice of a then for nonzero pressure that term 
leads to the tachyonic contribution to the gradient energy. However, after taking a 
closer look, it can be shown that the healthy contribution to the gradient energy, 
coming from the quartic Galileon, overwhelms the unhealthy one; rendering the 
background stable. A simple estimate shows that this remains to be the case as 
long as, p ^ p. For p ~ p, on the other hand, the coefficient of the gradient 
term of the perturbations is proportional to (p — 2p)/p, and for p > p/2 flips the 
sign to yield gradient instabilities. Such instabilities should manifest themselves in 
production of the tt quanta (albeit the coupling of these quanta to the matter is very 
weak) . Whether this may be dangerous, or interesting, for high density astrophysical 
objects, such as neutron stars, should be investigated including the fact that they're 
rotating and are neither exactly spherically symmetric, nor homogeneous; that is, 
they lack properties used in the derivation of above-mentioned results. 
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